arXiv: 1505.01580vl [physics.soc-ph] 7 May 2015 


Synchronization in an evolving network 

R K Singh^’0 and Trilochan Bagarti^’0 

^ The Institute of Mathematical Sciences, CIT Campus, Taramani, Chennai 600113, India 
^Harish-Chandra Research Institute, Chhatnag Raod, Jhunsi, Allahabad 211019, India 

In this work we study the dynamics of Kuramoto oscillators on a stochastically evolving network 
whose evolution is governed by the phases of the individual oscillators and degree distribution. 
Synchronization is achieved after a threshold connection density is reached. This cumulative effect 
of topology and dynamics has many real-world implications, where synchronization in a system 
emerges as a collective property of its components in a self-organizing manner. The synchronous 
state remains stable as long as the connection density remains above the threshold value, with 
additional links providing resilience against network fluctuations. 


Synchronization is a ubiquitous phenomenon observed 
across many natural and artificial systems. Flashing of 
fireflies, chirping of crickets, neurons in brain and Joseph- 
son junction arrays are well known examples that show 
synchronization. Kuramoto model provides a paradigm 
of ^nchronization of globally coupled phase oscillators 
[ij-lsj. However, in a system composed of a large number 
of components, no two of which have identical dynamics, 
synchronization emerges collectively in a self-organizing 
manner from a complex web of interactions amongst in¬ 
dividual components Q . Dynamics on comjdex networks 
provides a manifestation for such systems (g- In recent 
years a number of authors have studied the Kuramoto 
model on complex networks @ . Past studies have focused 
on static networks Q but a surge of interest in dynamic 
networks Q has led people to study the emergence of 
synchronization in dynamic networks In these stud¬ 
ies, the dynamics of the individual oscillators and the 
network topology has been treated independently. Syn¬ 
chronization emerges autonomously in many real systems 
which are constantly evolving in time. In this paper, we 
address the basic question of emergence of synchroniza¬ 
tion by providing a mechanism for the stochastic evolu¬ 
tion of the network governed by the system dynamics. 
Starting initially with a collection of independent oscil¬ 
lators we numerically show emergence of global synchro¬ 
nization from an interplay of network topology and local 
dynamics. We also find the synchronous state to be sta¬ 
ble as long as the number of nearest neighbors remains 
above a threshold with the additional links providing re¬ 
silience against network fluctuations M- 

We consider a system of N coupled phase oscillators 
described by the following equations : 

dO ^ 

^ = Wi -b — ^ ffy sin(6»j - 9i), i = (1) 

(Xl K/j, 

where 0i{t) and u>i denote the phase and frequency of the 
oscillator respectively, and e is the coupling strength. 
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Coupling between oscillators is described by the adja¬ 
cency matrix, g = [gij]. For the oscillator, = 
max{l,A:^}, where is the in-degree of the 

node i that counts the number of oscillators coupled to 
the oscillator. Similarly, out-degree, k~ = gji is 
the number of oscillators to which the oscillator is con¬ 

nected. The initial conditions are given by 0i(O) = 
for alH = 1,..., iV and g(0) = go. 

Evolution of the adjacency matrix g is governed by the 
dynamical rule : 

gij(tn+i) = 1 - gijitn) with probability P, 

tn-\-\ — tji -\- Tji, n — 0 , 1 , 2 ,..., ( 2 ) 

where P depends on k^, k^, 6i, 6j for i ^ j, and 
Tn are continuous random variables. From Eq. @ we ob¬ 
serve that the event of a single-flip in adjacency matrix 
is independent of previous flips (where a flip is defined as 
a change in value of gij from 0 to 1 and vice-versa). We 
make a further simplification by assuming that the flips 
occur at a constant rate. Hence, we choose to be expo¬ 
nentially distributed waiting time intervals of parameter 
A. Dynamics of the system proceeds in the following two 
distinct steps: (a) in a time interval oscilla¬ 

tors evolve according to Eq. o with adjacency matrix 
g(t) = g(tn) and (b) at time t = tn+i a pair of oscillators 
{i,j) chosen randomly are coupled or decoupled instan¬ 
taneously by the dynamical rule in Eq. Q- Note that 
due to stochastic evolution of the network by Eq. m, 
gij need not be symmetric. The mean waiting time 1/A 
and the average time period 2tt/ (w) provide two natural 
timescales of the system, where (uj) is the average of the 
frequency distribution. 

Two special cases of Eq. m and m are : (a) For 
gij = 1 forall i ^ j, and 0 otherwise, and A = 0 we obtain 
the original Kuramoto model [I|. (b) When the mean 
waiting time, A is finite and the transition probability P 
is constant, the steady state network is the well known 
Erdos-Renyi model pTl |. 

For the general problem defined by Eq. 0 and ([2]), the 

transition probability P for a given pair of oscillators i,j 

depends on the phases of 9i , 0, and the network degrees 
_|_ _ 1 _ 

of freedom, gij,kj^ , kj . We assume that the network de¬ 
grees of freedom are separable from the dynamical vari¬ 
ables 9i and 9j. The transition probability is dehned as 
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P = 9ij)VP +9^J{^P), ( 3 ) 

where 77 = r]{kf, k^) and p = /c(0i, 9j). Here, p is an even 
function of the phase difference 9i — 9j with a minimum 
at TT. We choose p as : 

p{9i,ej) = (2/7r)cos^((6»i -0j)/2). (4) 

We assume two different functions for 77 . Let 771 ( 2 ) = 1 
for kt = k~ = 0 and 

m {kt ,k~) = kj / (fc+ + kj ), (5a) 

m{kt, kJ) = l/{kf k~), (5b) 

for fc)*", kJ nonzero. Eq. (l5all is motivated by the obser¬ 
vation that hubs in a network are generally more influen¬ 
tial as compared to less interacting nodes Eq. (l5bl) 
on the other hand captures the form of dynamics where 
the tendency of a node to get connected to more nodes 
diminishes with its increasing neighbors [l3l| . 

Initially, all the oscillators are completely isolated (i.e. 
g(0) = 0) with the phases and frequencies distributed 
uniformly in the interval (0,27r) and (tt — 1/2, tt -|- 1/2) 
respectively. Average phase of the oscillators is defined 
as : 
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N ^ 

i=i 


idj 


( 6 ) 


where Oj is the phase associated with oscillator j and 7 / 
the average global phase. The magnitude r is a measure 
of the global synchronization of the system. Eor a com¬ 
pletely phase-locked state, r = 1. The connection density 
of the network is defined as : 

c =_ — _ (7) 

where E is the number of directed links in the graph at 
any given instant and N the total number of nodes, c 
is a measure of the denseness of the graph. Eor a com¬ 
pletely connected graph in which every node is connected 
to every other node, c = 1. 

The quantities of interest defined above, we proceed to 
solve the set of equations in Eq. (IT|) coupled with Eq. @ 
in order to study the evolution of the system as a whole. 
Numerical calculations suggest the existence of a critical 
value of global coupling strength, cq, below which the sys¬ 
tem remains in an asynchronous state but above which 
all the oscillators phase-lock with each other. We anal¬ 
yse the properties of the system for a value of coupling 
strength e > eo and then calculate the value of eo using 
finite size scaling M- 

Fig- Q] shows the variation of order-parameter, r and 
connection, c with time, t for global coupling strength, 
e = 1, for a system of size N = 100 with the parameter of 
the waiting time distribution, A = 1. As observed from 





FIG. 1: (Color online) Variation of r and c with t for the 
given initial conditions with e = 1, N = 100, A = 1 : (a) The 
system gets synchronized in almost the same time for both 
the transition rules Eq. (I5al) and Eq. (I5bl) . The inset shows 
the rise in order parameter for a static network topology at 
c = 0.02, 0.04. (6) c grows sublinearly with t, with the rate of 
increment different for both the rules, (c) Onset of synchro¬ 
nization in the system and its persistence after the network is 
not allowed to change any further once the system achieves a 
global phase-locking. If the links are removed randomly from 
the network then the synchronization is lost, (d) Connection 
density for the behavior in (c). The results are obtained after 
averaging over 10 ensembles. 


Fig. UKa), (inset), it takes a comparatively longer time for 
the system of phase oscillators to reach a phase-locked 
state in an evolving network as compared to a network 
of static topology with a given connection density. Sub- 
linear growth of c with t is shown in Fig. [Ijb). We have 
evolved the system upto t = 50000 and observed that the 
connection density c ~ t“,a < 1 for both the transition 
rules. The exact value of exponent a depends on the tran¬ 
sition rules and the value of coupling strength. The time 
series for r in the two figures also suggests towards exis¬ 
tence of a threshold value of connection density in order 
to achieve global synchronization. To confirm the exis¬ 
tence of a threshold value of connection density, we stop 
adding new links to the network once the system achieves 
global synchronization. We observe that the synchroniza¬ 
tion persists thereafter providing a strong indication to¬ 
wards the existence of a minimnm number of oscillators 
which should interact so that the system evolves coher¬ 
ently as a whole. Once the network stops evolving, we 
check the robustness of the synchronized state and find 
that synchronization is not resilient to random removal 
of links as the network gets fragmented(see Fig.[T](c),(d)). 

In Fig. [ 2 ] is the phase diagram r vs c show a mono¬ 
tonic increase of c with r. We observe that a threshold 
value of c is required in order to achieve synchronization. 
Once the synchronization is achieved, r becomes constant 
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FIG. 2: (Color online) Phase diagram of connection, c vs 
order-parameter, r. Once the synchronization is achieved, c 
rises abrnptly due to rise in network density, but more ad¬ 
dition of links has a negligible effect on the value of order 
parameter. 


but c rises due to addition of links on an average. Exis¬ 
tence of a threshold value of c can be understood from 
the observation that size of the largest connected compo¬ 
nent of the directed graph quickly becomes TV, the system 
size. This is achieved well before synchronization sets in. 
This also explains the desynchronization in the system 
due to random removal of links which makes the net¬ 
work fragmented into many components. This observa¬ 
tion was confirmed by walking on the graph using depth 
first search algorithm, which provides the size of the con¬ 
nected components in a graph M- A threshold value of 
c required for the system to achieve synchronization also 
suggests that a system of larger size will take longer time 
to achieve synchronization in comparison to a system of 
smaller size which is confirmed numerically. Numerically, 
exact value of threshold connection depends on the value 
of the order-parameter at which the system is considered 
to be synchronized, e.g., if r r; 0.9 is chosen to be the 
value of order-parameter for global synchronization, the 
corresponding threshold connection density is, c ~ 0.03. 

Fig. [3] shows the evolution of degree distribution with 
time for the network. The mean of the degree distribu¬ 
tion (fc^) grows sublinearly in time with {k~^) = {k~) 
for both the transition rules. The equality of mean de¬ 
gree is evident from the fact that total in-degree equals 
total out-degree for a directed network and hence, every 
node is affected almost equally by its neighbors. Vari¬ 
ances on the other hand are not identical for both the 
in- and out-degree distributions for rule 1, with the in¬ 
degree distribution showing less width in comparison to 
the out-degree distribution(Fig. [3](a),(b)). Note that for 
rule 1, rji{kf ,k~) ^ r]i{k~,k^), and since ? 7 i depends 
on the ratio of the out-degree of node j to the sum of in¬ 
degree of node i and out-degree of node j. Asa result dif¬ 
ferent oscillators contribute differently towards affecting 
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FIG. 3: (Color online) Degree distribution for the evolving 
network with N = 100, e = 1 and A = 1. (a) and (b) are the 
in-degree and out-degree distribution for the first transition 
rule of network evolution. Likewise (c) and (d) are for the 
second transition rule. The system was evolved upto t = 
50000 and the distribution curves are from t = 5000 to t = 
50000 in steps of At = 5000 with the curves shifting towards 
right in time. The results are obtained after averaging over 
10 ensembles. 
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FIG. 4: (Color online) Finite size scaling at mean-waiting 
time, A = 1 of r vs e curve for different system sizes, N for 
the two transition rules. The values of the critical exponents 
a and b are : a«0.5, 6«0.05. The critical value global cou¬ 
pling strength required to set the system of phase oscillators 
in a synchronous state is eo « 0.73. The results are obtained 
once the system reaches steady state. The calculations are 
done in t ^ oo limit and averaged over 10 ensembles. 


the dynamics of a given oscillator. However, for the sec¬ 
ond transition rule, we see that 'ri 2 {k^,k~) = 'r] 2 {k~,k^), 
therefore the width is identical for both the in- and out- 
degree distributions for obvious reasons. 

Finite size scaling of the transition curve(Fig.[3]) shows 
that the properties observed from numerical solution of 
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Eq. (U) do not depend on system size. Finite size scaling 
analysis also provides an approximate value of the crit¬ 
ical coupling strength, eq below which the system will 
never reach a phase locked state independent of its evo¬ 
lution time. The value of critical coupling strength is 
the same for both the transition rules and this can be 
understood in terms of the connection threshold. This 
is because once the system reaches synchronization, in¬ 
creasing connection density further does not have any 
significant effect on the order parameter. Furthermore, 
a finite value of critical coupling strength is also evident 
from the observation that the second moment of the de¬ 
gree distribution of the steady-state network is finite [l^ . 

We benchmark our results against existing literature 
by solving the dynamical equations with oscillator fre¬ 
quencies distributed according to Cauchy distribution 
peaked at the origin and having scale parameter, 7 = 
0.25. In the mean field limit for a completely connected 
graph, the critical value of coupling strength is, eg = 27 
[3|. We evolve the network by adding links with proba¬ 
bility one to randomly chosen pair of nodes so that the 
steady state topology of the network is that of a com¬ 
pletely connected graph. We find eg ss 0.4, which pro¬ 
vides a benchmark to our results. 

In order to investigate the dependence of eo on the 
mean of the waiting time distribution A, we looked at the 


variation of order-parameter with coupling for different 
values of A for a fixed system size. We find that varying 
A does not have any appreciable effect on the value of 
critical coupling strength for either of the rules. This 
observation suggests that eg is independent of A. 

To conclude, we have shown that the interdependence 
of topology and dynamics is necessary in order to un¬ 
derstand how synchronization emerges collectively from 
local interactions in a self-organizing manner. This has 
real world implications for systems in which interactions 
keep changing with time but an overall rythm persists on 
an average. It is also observed in natural systems that the 
synchronous state does not persist forever. In this model 
we have provided a mechanism for collective emergence 
of synchronization from local dynamics. The synchro¬ 
nization achieved is characteristic of the network evolved 
by the transition rules. However, it is independent of the 
connection density of the network provided it is above 
a threshold. Existence of a connection density threshold 
in order to achieve synchronization has practical applica¬ 
tions, e.g., from the context of Josephson junction arrays 
in reducing the number of redundant connection links in 
order to achieve desired output thereby minimizing the 
cost of wirings. We hope that this work provides an im¬ 
petus for future work on synchronization from a dynamic 
perspective. 
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